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ABSTRACT 

We compute the profiles of resonance doublet lines (S'1/2 — ^*1/2,3/2) formed in bipolar 
winds with velocity greater than the doublet separation in symbiotic stars. Particular 
attention has been paid on the doublet line ratio, where an essential role is played 
by the conversion of the short wavelength component arising from the S'1/2 ~ A/2 
transition into the long wavelength component for the transition S1/2 — A/2- We 
adopted a Monte Carlo technique and the Sobolev approximation. Our bipolar winds 
take the form of a cone and are characterized by the terminal wind velocity, the mass 
loss rate and the opening angle of the cone. When an observer is in the polar direction 
and the Sobolev optical depth rsob — 1 , we mainly obtain profiles with inverted flux 
line ratios, where the short wavelength component is weaker than the long wavelength 
component. When an observer is in the equatorial direction, we find that the prohles 
are characterized by two broad components, where the long wavelength component is 
the broader and stronger of the two. We conclude that the prohles obtained in our 
model provide a qualitative understanding of broad prohles and inverted intensity 
ratios of the doublets in symbiotic stars. 

Key words: line: formation — line: prohles — radiative transfer — scattering — 
binaries: symbiotic — stars: winds, outhows 


1 INTRODUCTION 

Symbiotic stars are generally known to be interacting bina¬ 
ries consisting of a red giant (or a Mira-type variable) and 
a hot companion that is usually a white dwarf (e.g. Kenyon 
1986). Most red giant components suffer a heavy mass loss in 
the form of a slow stellar wind with a typical terminal speed 
10 ~ 20 km s“^ that is comparable to the escape velocity of 
a giant (e.g. Schmid 1996). In contrast, many white dwarf 
systems including planetary nebulae are known to possess 
fast outflows with terminal speed > 1000 km s“^. Similar 
fast winds are also known in some symbiotic stars including 
AG Peg (Vogel & Nussbaumer 1994). 

The orbital elements of symbiotic stars are not well- 
constrained, but light curves often show that they possess 
a long period of several hundred days (e.g. Iben & Tutukov 
1996). Miirset & Schmid (1999) presented the relation be¬ 
tween the spectral types of the cool giants and the orbital 
periods to reveal that almost all symbiotic stars are well- 
detached binary systems. Therefore the most important bi¬ 
nary activity may be found from the interaction between 
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the two different kind of winds. These two winds may col¬ 
lide forming a shocked region, which can be identified with 
X-ray observations and is consistent with the fact that sev¬ 
eral symbiotic stars are known to be X-ray sources (Girard 
& Wilson 1987, Miirset, Wolff & Jordan 1997, Ezuka, Ishida 
& Makino 1998) 


In many symbiotic stars the Raman-scattered O vi 
AA6827,7088 features exhibit multiple peak profiles and 
strong polarization accompanied by the polarization flip 
(Schmid 1989, Schmid & Schild 1994, Harries & Howarth 
1996). Lee & Park (1999) proposed that these features can 
be explained by assuming that there is an accretion disk 
around the white dwarf formed through capture of the slow 
wind from the giant (e.g. Mastrodemos & Morris 1998). Ac¬ 
cording to the theoretical modelling by Paczynski & Zytkow 
(1978), periodic hydrogen shell flashes may occur in a white 
dwarf with the accretion rate 10“^^ — 10“^ Mq yr“^, where 
each outburst may last for decades. With this eruption the 
radiative pressure will drive a stellar wind around a white 
dwarf, which may take a bipol ar form subject to the cir - 
cumstellar matter distribution ( Soker fc Rapparport 2000 ). 
This interpretation is interesting because most sym biotic 
stars with known nebular morphologies are bipolar (Cor- 
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radi 1995). Currently it is very controversial whether bipolar 
planetary nebulae possess central binary systems (e.g. Soker 
1998). 

It is, therefore, very interesting and important to in¬ 
vestigate the line profiles that indicate the fast outflowing 
motion, from which we may find the physical properties 
associated with the bipolarity of the wind. Since the res¬ 
onance doublets arise from the common electronic transi¬ 
tions S'i /2 — Pi/ 2 , 3/2 and have the separation ranging from 
500 km (for C iv) to 1,650 km s“^ (for O vi), these 
lines can be an excellent tool to investigate the outflowing 
hot wind around the white dwarf component. 

When the bulk velocity changes significantly in a re¬ 
gion that is much smaller than the scale height of the phys¬ 
ical quantities such as density, the radiative transfer can be 
described by the Sobolev approximation. In this case, the 
optical depth for a line photon is inversely proportional to 
the velocity gradient in the direction of the photon prop¬ 
agation (e.g. Sobolev 1947; Rybicki & Hummer 1978). In 
the case of resonance doublet lines, line photons arising 
from the 51/2 — P3/2 transition will be resonantly scattered 
with 51/2 — Pi/2 transition by receding ions with the speed 
of the doublet separation, if the outflow is an accelerating 
wind with a speed larger than the doublet separation. We 
denote this type of scattering by ‘double scattering.’ This 
double scattering converts line photons associated with the 
S'1/2 — P3/2 transition into line photons associated with the 
S'1/2 — Pi /2 transition. This will change the intrinsic doublet 
line ratio, which is 2 : 1 in the optically thin limit and 1 : 1 
in the optically thick limit. Olson (1982) investigated this 
problem and applied to stellar winds around O and B stars. 
His main concern was limited to the investigations of various 
P Cygni profiles in spherical winds. 

lUE observations indicate that many objects including 
symbiotic stars and planetary nebulae show various doublet 
line ratios between 2 : 1 and 1 : 1 (Feibelman 1983, Schmid 
et al. 1999). The microphysical processes that lead to these 
various line ratios may be associated with the existence of 
dust and/or collisional de-excitation (Ahn & Lee 2002, in 
preparation). 

Michalitsianos et al. (1988) showed that some symbi¬ 
otic stars including RX Pup and R Aqr show anomalous 
line ratios, in which the short wavelength component of the 
C IV doublet was observed to be weaker than the long wave¬ 
length component. Vogel & Nussbaumer (1994) showed that 
the symbiotic nova AG Peg exhibits broad He ii emission 
lines that are formed in a fast wind. In this system, the res¬ 
onance doublets N v and C iv exhibit inverted line ratios 
where the short wavelength component is weaker. The au¬ 
thors proposed that the line formation is strongly affected 
by the double scatterings, whereby signihcant fraction of 
photons are converted. Furthermore, the doublets did not 
display P Cygni profiles. These facts imply that the out¬ 
flowing motion is confined to specific directions (plausibly 
in the polar directions), which exclude the observer’s line of 
sight. 

There have been many theoretical investigations on the 
P Cygni profiles shown in the resonance lines of metal el¬ 
ements by adopting the Sobolev approximation. However, 
relatively little attention has been paid on the profiles ob¬ 
served outside the wind flowing direction. In this paper, we 
perform Sobolev Monte Carlo computations to obtain the 


profiles of resonance doublet lines formed in bipolar winds 
that may be present in symbiotic stars. 

In Section 2, we briefly describe the Sobolev theory 
and the basic atomic physics concerning the resonant dou¬ 
blet lines. We also present the kinematic stellar wind model 
adopted in this work and the Monte Carlo procedure. Our 
results are presented in Section 3. Finally, we summarise and 
discuss our results and observational implications in Sec¬ 
tion 4. 


2 WIND MODEL AND CALCULATIONS 


2.1 The Sobolev Theory 


In the case of the stellar wind around the white dwarf com¬ 
ponent in a symbiotic star, the wind is highly ionised and 
characterized by the thermal speed ^ 10 km s“^ and the 
bulk speed of order 10® km s“^ (Girard & Wilson 1987). 
We may apply the Sobolev approximation to describe the 
radiative transfer, if the acceleration of the wind material 
occurs in a small region where there are no considerable 
changes in physical quantities. We will check the validity of 
the Sobolev approximation later in detail. 

When a scattering medium is of moderate column den¬ 
sity, wing scatterings can be safely neglected. In this case, 
the scattering cross section associated with a resonance line 
photon can be well-approximated by the Dirac d-function 
with the strength fabsi'^e'^/mec), or 


— fabs 5{y r'o), 

meC 


( 1 ) 


where fabs is the absorption oscillator strength associated 
with the resonance transition and the frequency ly is mea¬ 
sured in the rest frame of the scattering atom. Therefore, a 
scattering occurs only when the resonance condition is met. 
Because of the local thermal motion, the physical width of 
the scattering region of a photon is given by the width where 
the difference of velocity is of order the thermal velocity Vth, 
that is, Ar ~ vth{dv/dr)~^. 

Noting that the local velocity distribution of scattering 
atoms is given by the Maxwell-Boltzmann distribution, a 
straightforward computation gives the Sobolev optical depth 
Tsob for a given direction. 


'T~Sob — 


advc/[v{dV/ds)] 


/^e^N 


<U 

\ds) 
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where s is the distance along the photon propagation direc¬ 
tion, Ao represents the resonance wavelength, and Ui is the 
ion number density. A more complete and detailed descrip¬ 
tion of the Sobolev theory can be found in the literature 
(e.g. Sobolev 1947, Rybicki & Hummer 1978, Lee & Bland- 
ford 1997, Ahn, Lee & Lee 2000). 

In this work, we are particularly interested in the reso¬ 
nance doublets, where the doublet separation is smaller than 
the wind velocity. Those doublets include C iv AA1548,1551, 
N v AA1239, 1243, and O vi AA1032,1038, which are known 
to be prominent in symbiotic stars from observations made 
with space instruments such as lUE, the Hopkins Ultra¬ 
violet Telescope (HUT), and the Orbiting and Retrievable 


© 0000 RAS, MNRAS 000, 000-000 



Resonance Doublets in Bipolar Winds 3 



1 0,8 0.6 0.4 O.Z 0 

1/x = R./r 


Figure 1. The velocity field normalized by the terminal velocity 
and the corresponding Sobolev optical depth rgoj, according to 
the distance normalized by i?*, the radius of the photosphere of 
the hot white dwarf. Solid and dashed lines represent the Sobolev 
optical depths of the short and the long wavelength components 
of the resonance doublets, respectively 


Far and Extreme Ultraviolet Spectrometer (ORFEUS). The 
Si /2 — P 3/2 transition corresponding to the short wavelength 
component has the twice stronger oscillator strength than 
the Si /2 — Pi /2 transition for the long wavelength compo¬ 
nent. Therefore, at the same position, the Sobolev optical 
depth Tg^j, corresponding to the short wavelength compo¬ 
nent is also twice the optical depth for the long wave¬ 
length component as are shown in Fig. From now on, we 
refer the Sobolev optical depth rsob to be Tg^;, when we do 
not specify the transition. 


2.2 Kinematics of the Fast Stellar Wind 

The existence of fast stellar winds with speed v ~ 
10® km s“® around a white dwarf has been known in many 
white dwarf systems including the central stars of plane¬ 
tary nebulae (e.g. Cerruti-Sola & Petrinotto 1989). Many 
P Cygni profiles formed in fast stellar winds have been suc¬ 
cessfully fitted using a velocity law v{r) = i!oo(l — r/Rt)^ 
with the Sobolev approximation (Tamers, Cerruti-Sola & 
Perinotto 1987), where the stellar wind is assumed to start 
at the photospheric radius i?*, /3 is a dimensionless param¬ 
eter and Voa is the terminal velocity. We scale the physical 
distance with i7* so that the distance r from the stellar cen¬ 
tre is obtained from the dimensionless parameter x defined 
by 

X — r/Rt. (3) 

The stellar wind in this work is assumed to be steady 
and have geometry in the form of a bipolar cone includ¬ 
ing the spherically symmetric case. The velocity v{r) of 
the wind increases monotonically outward and asymptoti¬ 
cally approaches the terminal speed Voo- Although the exact 


model of the velocity law in symbiotic stars is unknown, 
Castor, Abbott & Klein (1975) proposed a very steep ve¬ 
locity law for a radiation-driven stellar wind. We measure 
the wind velocity with the terminal speed, and introduce a 
dimensionless velocity w defined by 

■W = v/voo- (4) 


Instead of assuming a wind velocity law, we prescribe 
the Sobolev optical depth by 



where tq is the initial Sobolev optical depth at the wind base 
and e is a dimensionless positive number. We choose e to be 
small, so that the Sobolev optical depth does not decrease 
significantly in most of the wind region. This will enhance 
the effect of double scattering, where continuum photons 
blueward of the S 1/2 — P 3/2 transition may get scattered by 
the transition S 1/2 — P 1/2 further downstream of the wind. 
We use the initial Sobolev optical depth rg^;, q = 3, Q — 
1.5 for the short and the long wavelength components of the 
doublet. 

The density profile n(r) is obtained from the mass flux 
conservation assuming the constant mass loss rate M, 


n{r) 


M 

4:71 mpr^v{r)cf 


1 

r^v{r) ’ 


( 6 ) 


where mp is the proton mass and c/ is the covering factor 
of the wind normalized by the whole sky. From Eqs. (2), (5) 
and (6), the velocity law of the wind is given by 




-l-|-e-| 1/2 


( 7 ) 


In order to check the validity of the Sobolev approximation, 
we consider the length scale over which the bulk velocity 
changes by the amount of Vth, 


Lsob 


Vth 

dv I dr 


( 8 ) 
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'^00 



We note that Vthivoo < 10“^, and therefore Lgob < 
0.0277**^“". This implies that Lgob can be comparable to 
Rt, when a: > 10 and e = 0.3. When x is small, the den¬ 
sity changes steeply near 77*, and Lgob 77* validating the 
use of the Sobolev approximation. When x gets larger, the 
density variation is almost scale-free, which also makes the 
Sobolev approximation reliable. 

The mass loss rate is related to the Sobolev optical 
depth TSob by 


M = 5.1 X 10-Woil - e)cf Mo yr-\ (9) 


and the hydrogen number density can be obtained by 
n(r) = 1.6 X 10^ro(l — e)x~^w~^ 

^ (lOiicm) (l03 km s-i) 

According to the study of Cerruti-Sola & Perinotto 
(1989), the terminal velocities of the fast winds in many 
planetary nebulae range 1000 — 2500 km s“^. Vogel & Nuss- 
baumer (1994) also showed that the fast wind in the sym¬ 
biotic star AG Peg is about 1000 km s“^. In this study, 
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Figure 2. The spectrum of the source that consists of a flat con¬ 
tinuum and resonance doublet lines, which we assume is N v here. 
We used Gaussian profiles for the line doublet whose equivalent 
widths are 5 A and 2.5 A for the short and the long wavelength 
components, respectively. The velocity dispersion is chosen to be 
(T = 30 km s“^ for both components. The wavelength regions are 
divided according to the wavelength relative to the resonance line 
centres of the doublet. 

we choose the doublet separation to be 960 km s“^ corre¬ 
sponding to the N V doublet and set the terminal velocity 
Voo = 3000 km s“^ in order to investigate the effect of dou¬ 
ble scattering. Therefore, we focus on the prohles formed by 
various scatterings, not on the full dynamical range of the 
wind. 

With the choices of the parameters e = 0.3, tq = 3 
and the wind half opening angle 30°, the mass loss rate 
M = 1.3 X 10“ Mq yr“^. It is notable that the mass loss 
rate is compatible to those of Michalitsianos et al. (1988) 
which are 4.5x 10“^^ < M < l.Ox 10“^^ Mq yr“^ for R Aqr 
and 1.5 x lO”*^^ L M < 3.1 x lO”*^^ M© yr“^ for RX Pup. 
According to them, it is estimated that the momentum flux 
from the red giant is about 200 times greater than that from 
the hot companion. The corresponding mass loss rate from 
the red giant is 1.3 x 10“^ Mq yr“^ which is consistent with 
the calculation by Vogel & Nussbaumer (1994). 

The profiles of the wind velocity and the Sobolev optical 
depth are shown in Fig. as a function of the normalized 
distance x. 

2.3 The Sobolev Monte Carlo Code 

Firstly, we consider the volume emission measnre in the stel¬ 
lar wind region. The relevant quantity to be considered is 

J n{rfdV = nlRlTS{l-e)A7T 

= noR2To(l - e)^47r/(e), (11) 

where no = 1-6 x lO’’ cm“® and /(e) = 

which is a steep fnnction of e near zero. For a small value 



Figure 3. Conceptual diagram for single resonance scatterings 
and double resonance scatterings. The solid, dotted, dashed and 
wavy lines represent no scattering, single scattering resonant with 
Si /2 — ^ 3 / 2 ! single scattering resonant with S 1/2 ~ P 1/2 and dou¬ 
ble scattering, respectively. The indices represent possible scat¬ 
terings that a given photon in Region I can suffer. 


e = 0.3, we also have a moderate value of /(e) and the 
quantity n^V is basically determined by noR* — lO'*^ cm“^ 
much smaller than typically observed in symbiotic stars (e.g. 
Proga, Kenyon & Raymond 1998). Therefore, if we consider 
the wind as a photon scattering region, almost all the pho¬ 
tons are generated near the hot star component. This means 
that the photon source is separated from the wind. 

From this consideration, we prepare a photon source 
near the white dwarf component, from which both contin¬ 
uum photons and line photons are injected into the wind 
region. We choose the equivalent width of the short wave¬ 
length component to be 5 A and 2.5 A for the long wave¬ 
length component. This ratio corresponds to the optically 
thin limit of the emission plasma. Without specifying the 
details of the kinematics of the emission region, we just as¬ 
sume that the line proHle is given by a Gaussian with the 
velocity dispersion a = 30 km s“^ for both components. 
Fig. 1^ shows the proHle of the incident photon flux that is 
injected into the wind region. 

According to the wavelength relative to the two reso¬ 
nance line centres, we conceptually divide the injected pho¬ 
tons into three categories, which we call Region I, II and III 
shown in Fig. 

For photons in Region III, there is no point in the wind 
where the resonance condition can be satisfied for either 
transition. Therefore, no photons in Region III will be scat¬ 
tered in the wind in our model. Photons from Region II can 
only be resonantly scattered in the wind with the S1/2 — R1/2 
transition, for which the relevant Sobolev optical depth is 
0.5rso6- 

Most complicated interactions can be seen for photons 
in Region I, because they can be first resonantly scattered 
with the Si/2 —R3/2 transition and can be additionally scat- 
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X = r/R. 



X = r/R. 


Figure 4. Examples of the Sobolev surfaces for double scattering 
in the x — y plane. A given photon is resonantly scattered at the 
position marked by an asterisk, and the surfaces with the velocity 
component equal to the doublet separation. The arrows and the 
solid lines represent the velocity vectors at the positions and the 
wind region, respectively. 



6 = cos^/fi • f) 

Figure 5. The dependence of the Sobolev optical depth Ts^b 
on the photon propagation direction at a given position. The 
horizontal axis represents the angle between the incident pho¬ 
ton wavevector and the radial or the wind direction. The solid, 
dotted and dashed lines represent the inverse of the velocity gra¬ 
dient at x=1.16, 1.2 and 1.3 which are normalized by the inverse 
value of the velocity gradient for the radial direction. 


tered with the S1/2 — P1/2 transition. Hence there are four 
cases that can occur for photons in Region I. The first case 
denoted by I-l in Fig. ^ corresponds to the escape without 
any scattering. The second case (T2) represents the escape 
after a resonance scattering associated with S1/2 — P3/2 tran¬ 
sition. The third case (T3) corresponds to the escape after a 
resonance scattering associated with S1/2 — P1/2 transition. 
The final case (T4) is the double scattering, where a given 
photon escapes from the bipolar wind after resonance scat¬ 
terings associated with both transitions. Fig. ^ illustrates 
these four possible cases. 

We consider the Sobolev surface, which is defined as the 
collection of those points moving away from a given point 
with the same velocity component along the line of the pho¬ 
ton propagation. In order to describe doubly scattered pho¬ 
tons, we have to consider the Sobolev surface moving away 
with the velocity of the doublet separation from a given 
point where the first resonance scattering occurs. For illus¬ 
tration, in Fig. ^ we show a couple of Sobolev surfaces for 
double scattering in x — y plane, where the given photon 
is firstly resonantly scattered at the positions marked by an 
asterisk. The Sobolev optical depth associated with this 
double scattering is given by 


T^ob = nifLXo( — )\n ■ (n • V)yr 
xrrieCj 

= ni/abs^of^^) {n- 

\meCj 


■ -I- — [1 - (n • f) ] 
dr r 


( 12 ) 


for various directions at three positions x = 1.16,1.2, 1.3 in 
order to find the dependence of on the photon propa¬ 
gation direction. We can see that ma-y vary by a factor 
of three dependent on the direction. 

Since we are particularly interested in bipolar winds, 
which may be related with the bipolar nebular morpholo¬ 
gies known in a number of symbiotic stars (e.g. Corradi 1995, 
Corradi & Schwarz 1993), it is natural to expect that there 
exists a very thick circumstellar component in the equato¬ 
rial plane. Therefore, in our model, we introduce an opaque 
medium outside the photon source and the wind base. For 
simplicity, we set the physical size of the opaque medium to 
be X = 1.01, so that observer in the polar direction can see 
the wind base, which is hidden from the equatorial direction. 

We typically inject 10® continuum photons in a bin with 
A A = 0.1 A. The injection number is adjusted near the 
doublet line centres so that the line features possess the 
equivalent widths of 5 A and 2.5 A and the Gaussian profiles 
with velocity dispersion u = 30 km s“^. We collect emergent 
photons according to the direction cosine of the wave vector 
with a bin size of A/i = 0.1, and the wavelength bin AA = 
0.023 A corresponding to the thermal width of T = 10^ K. 
In our simulations, we normalize the relevant quantities to 
those associated with N v AA1238.8, 1242.8. 


3 RESULTS 


where n is the photon wave vector before double scatter¬ 
ing and f^bs is the oscillator strength corresponding to the 
Si/2 — Pi/2 transition. 

In Fig. 1^ we show the inverse of the velocity gradient 


We investigate the dependence of the profiles on the ob¬ 
server’s line of sight, the mass loss rate, the wind half open¬ 
ing angle and the existence of the opaque circumstellar re¬ 
gion. 
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Figure 6. Profiles according to the various observer’s lines of sight. The number of photons counted in each observer’s line of sight is 
referred to the relative flux and is shown as a solid line. The dotted and dashed lines mean the number of photons which suffered single 
scattering and the double scattering, respectively. The wind half opening angle is set to be 45°. 


3.1 Profiles for Various Lines of Sight 

We first fix the wind half opening angle to be 45° and con¬ 
sider the variation of the profiles dependent on the observer’s 
line of sight. Fig. ^ shows our result, in which we find signifi¬ 
cant difference in the profiles viewed from the polar direction 
and from the equatorial direction. 

When TSob is much larger than 1 and the observer is in 
the polar direction, we normally obtain P Cygni type pro¬ 
files. However, when the Sobolev optical depth rsob is mod¬ 
erate {rsob — 1), the observer sees profiles characterized by 
two narrow peaks. The blue part of the short wavelength 
peak is obtained from suppression of the incident flux by 
a factor e“'’'so!>(''i)“0.5-rg„6(r2) .^^^ith ri and r 2 being the lo¬ 
cations of the resonance transitions with S 1/2 — P 3/2 and 
Si /2 — Pi/ 2 , respectively. In contrast, the red part of the 
short wavelength peak and the blue part of the long wave¬ 
length peak are suppressed by a factor where rs 

is the location of the resonance scattering with S 1/2 — Pi/ 2 - 
Therefore the double scattering process affects the doublet 
line ratio significantly. In our model, only the blue part of the 
long wavelength component can satisfy the resonance con¬ 
dition and be affected by the wind. Hence we observe the 


fairly stronger long wavelength component than the short 
wavelength. This is interesting considering the lUE spectra 
of RX Pup and R Aqr which were pointed out by Michal- 
itsianos et al. (1988). According to them, in RX PUP and 
R Aqr the C iv doublet showed anomalous line ratio in the 
blue part only, where the short wavelength component is 
weaker than the long wavelength component. However, in 
the red part, the doublet flux ratio almost restores the op¬ 
tically thick limit of 1:1. Therefore, it is very probable that 
the doublet line ratio in the blue part may be affected by 
extrinsic components including the fast outflowing wind. 

As is apparent from Fig. the continuum blueward 
of the short wavelength component (Region I in Fig. ^ is 
weaker than that for Region H, which is again weaker than 
that for Region HI. A similar trend is observed for fj, > 0.85 
(in the direction of the wind region). However, qualitatively 
different profiles are obtained for p < 0.65, for which the 
observer is outside the wind flowing direction. 

If we see symbiotic stars outside the wind flowing direc¬ 
tion (/i < 0.65 in Fig. ^), the emergent profiles are composed 
of only scattered photons and no transmitted continuum can 
be seen. One of the most interesting features of this bipo¬ 
lar wind system is shown for /r ~ 0, where we see the clear 
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Figure 7. Profiles for various mass loss rates M. We fix the observer’s line of sight to be in the polar direction. The profiles shown in 
the left correspond to spherical winds. In the right part are shown the profiles for bipolar winds with the half opening angle 30°. The 
mass loss rates are chosen so that the initial Sobolev optical depth tq = 0.6,3 and 15 and Cf represents the covering factor of the wind 
with respect to the photon source. The various lines represent the same quantities illustrated in Fig. 0 


effect of the photon redistribution by the double scattering 
process. The net effect of double scattering is to take out the 
flux of the short wavelength component and redistribute it 
around the long wavelength component. Therefore, the flux 
around the long wavelength component increases. 

More quantitatively, the ratio of the fluxes contributed 
by single scatterings at the short and the long wavelength 
components and double scattering to the total flux is given 
by 15%, 39% and 46% in the particular case 0.15 < /r < 0.25 
in Fig. This means that more than half of 1239 pho¬ 
tons are doubly scattered and converted into 1243 photons. 
It is notable that the profile around the long wavelength 
component gets greater dispersion than the dispersion that 
would be obtained in the case when only single scattering 
at the long wavelength component is operational. This is 
because doubly scattered photons acquire the Doppler shift 
corresponding to the doublet separation in addition to the 
Doppler shift for the velocity at the scattering point. There¬ 
fore, the doublets in a bipolar wind are characterized by the 
inverted flux ratio and much broadened profiles of the long 
wavelength component. 


3.2 Dependence on the Mass Loss Rate 

In order to investigate the dependence of the profiles on the 
mass loss rate M, we set the covering factor of the wind 
to be 1 and 0.13 corresponding to the spherical and bipolar 
winds with half opening angles 90° and 30°, respectively. In 
this subsection, we fix the observer’s line of sight to be the 
polar direction. 

According to the various mass loss rates, considerably 
different profiles are obtained as shown in Fig. |^. The mass 
loss rate is proportional to the initial optical depth to as in 
Eq. (9). Various doublet line ratios are obtained because the 
peaks of the resonance doublet are sensitive to tq. For the 
spherical wind, the Sobolev surface for double scattering at a 
given position subtends a fairly large solid angle, and there¬ 
fore the emergent flux profile will be significantly affected 
by photons redistributed by double scattering. In particu¬ 
lar, when the mass loss rate is very high (the top-left panel 
of Fig. 1^, the suppression of the short wavelength compo¬ 
nent is so large that no apparent line feature is found around 
at 1238.8 A. Furthermore, all the photons blueward of the 
S 1 / 2 —P 3/2 resonance are redistributed in frequency space up 
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Figure 8. The profiles for the wind with half opening angle 60° and 30°. The various lines represent the same quantities illustrated in 
Fig.§ 


to the wind terminal speed and form the strong and broad 
component around the S1/2 — P1/2 resonance wavelength. 

However, for the narrow bipolar winds, the Sobolev sur¬ 
face for double scattering has a small solid angle, where we 
expect a negligible contribution from doubly scattered pho¬ 
tons. Since the scattering process in the bipolar wind can be 
effectively regarded as a photon annihilation process from 
the observer’s line of sight, we see a quick disappearance of 
the broad components around the doublet formed by the 
photons redistributed by single scattering. 

3.3 Dependence on the Wind Opening Angle 

We briefly investigate the profile dependence on the wind 
opening angle in order to study the profile formation in a 
bipolar wind. As a reference, the profile from a spherical 
wind can be found in Fig.^. The inverted intensity ratio 
will be sensitively dependent on the mass loss rate. Around 
the line centres, there appear broad bases which are formed 
by those photons redistributed from single resonant scat¬ 
terings associated with S1/2 — Pi/2,3/2 transitions, and also 
contributed from double resonant scatterings of photons ini¬ 
tially generated in Region I. 

In Figs. ^ we generated profiles for bipolar winds with 


half opening angles 60° and 30°, respectively. We chose 
two lines of sight for illustration. The overall trend is ex¬ 
plained in the previous subsection. As the half opening an¬ 
gle decreases, the redistributed photons contribute less to 
the profile formation. This effect is clearly seen in the case 
0.25 < fj, < 0.35, for example. In a bipolar wind with 60° half 
opening angle, we can barely recognise the short wavelength 
component, whereas the long counterpart is enhanced signif¬ 
icantly. In this case, there is effectively one broad line com¬ 
ponent that is centred nearly at 1242.8 A, slightly blueward 
of the long wavelength component. However, in the case of 
the 30° bipolar wind, there appear two clear line components 
of which the longer one is the stronger. The contribution of 
redistributed photons is significantly suppressed. The long 
wavelength component is, however, broader than the short 
wavelength component and there is also a broad base red- 
ward of the long wavelength component. 

The fluxes observed from the equatorial direction dis¬ 
play various profiles that are not of a P Cygni type. These 
various profiles are typified by inverted line ratios, and are 
formed by photons redistributed by single and/or double 
scatterings, which play an essential role in profile formation. 
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Figure 9. Profiles for the cases where the opaque media are changed. In the left, the opaque media are shrunken to hide only the photon 
injecting source and the wind base is fully seen from all directions. The opaque media are completely removed in the right, so that the 
photon injecting source and the wind base are fully seen. The wind half opening angle is set to be 45°. The various lines represent the 
same quantities illustrated in Fig. 0 


3.4 Miscellaneous Cases 

In Fig. we shrink the opaque component, so that it only 
hides the photon source. By doing this, effectively we con¬ 
sider the case where the wind base part is fully seen from 
any directions. 

When an observer is in the polar direction, we get the 
same result as before. However, from the equatorial direc¬ 
tion we obtain different profiles to which the wind base part 
contributes additionally with very small Doppler shifts. For 
the equatorial viewer, the wind base part actually acts as 
a narrow line source which is almost at rest. Therefore, the 
qualitative difference is due to the addition of the narrow 
line component to the broad scattered profiles. 

Finally, when we remove the opaque component 
(Fig. then we see fully the source and the scattering 
wind. No discernible difference is found for profiles observed 
from the polar direction. For an observer in the equatorial 
direction, the profile is dominantly that of the source and 
modified in a very small amount by scattered components. 


4 DISCUSSION AND SUMMARY 

In this paper we investigate the doublet resonance line for¬ 
mation in bipolar conic flows in symbiotic stars by adopting 
the Sobolev Monte Carlo method. We also pay special at¬ 
tention to the doublet line ratios. 

When an observer is in the polar direction, there ap¬ 
pear two peaks at the line centres coinciding with those of 
the resonance doublets. This is because we chose the mass 
loss rate M so that the Sobolev optical depth rsob — 1 for 
most of the velocity space, and therefore the line flux of the 
short component is suppressed by a factor ~ 0.3. 

The red part of the long wavelength peak does not suffer 
any suppression. This fact implies that the detailed peak 
profiles differ in a systematic way. The suppression pattern 
also repeats for continuum parts, where the weakest occurs 
blueward of the short wavelength peak. 

Michalitsianos et al. (1988) investigated the profiles of 
C IV AA1548,1551 in the symbiotic stars RX Pup and R Aqr, 
from which they found the weaker short wavelength compo¬ 
nent than the long wavelength component. They particularly 
noted that the profiles are characterized by significant sup- 
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pression in the blue part of the short wavelength peak. The 
peak flux ratio is most sensitively dependent on the Sobolev 
optical depth, particularly when rsob — 1 and is insensitive 
on the covering factor of the wind. 

In RX Pup and R Aqr, there may exist a scattering 
component outflowing with speed less than the doublet sep¬ 
aration. In this case, only single scattering is operational, 
which can cause the deviations of the doublet line ratio from 
the optically thin limit. However, when double scattering is 
possible, much more various doublet line ratios can be pro¬ 
duced. It is uncertain that there exists an outflow moving 
faster than the doublet separation of 500 km s“^ in RX Pup 
and R Aqr. In order to find the concrete evidence of double 
scattering, we need to obtain spectra with sufficient quality 
to discern the continuum level which are contributed by the 
redistributed photons. 

When we observe in the equatorial direction, we find 
that the profiles are characterized by two broad components, 
where the long wavelength component is the broader and 
stronger of the two. This qualitative feature does not change 
even if we vary either the mass loss rate or the Sobolev 
optical depth Tsob as long as Tsob exceeds unity. However, 
when the observer is in the polar direction, various P Cygni 
type profiles are obtained dependent on the mass loss rate. 
We believe that the bipolarity of the stellar wind is mainly 
responsible for these spectroscopic characteristics. 

It is quite interesting that the symbiotic nova AG Peg 
shows the profiles of N v and C iv in the hot wind phases 
similar to those obtained in this work. We strongly believe 
that the hot wind in AG Peg takes the bipolar form with 
the speed exceeding 10® km s“^. Kenny, Taylor & Seaquist 
(1991) used the Very Large Array (VLA) in order to inves¬ 
tigate the morphological structure of AG Peg. They found 
a bipolar outer nebula extending ^ 40 They also found 
a bipolar enhancement in the inner nebula having a subarc- 
second extent, which is consistent with the existence of a 
bipolar wind. However, it should also be noted that Nuss- 
baumer, Schmutz & Vogel (1995) reported a P Gygni type 
absorption in N v in AG Peg using the Hubble Space Tele¬ 
scope. 

Most symbiotic stars are not optically resolved and 
therefore their morphology is very uncertain. In planetary 
nebulae, it is still highly controversial whether the binarity of 
the central star system can be linked to the bipolarity of the 
nebula (e.g. Morris 1987, Livio & Soker 1988, Soker 1998). 
A significant input should be provided from the morpholog¬ 
ical study of symbiotic stars. In this respect, it is proposed 
that the resonance doublets can be a useful diagnostic of the 
bipolar winds in symbiotic stars. 
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